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Abstract 

Due to the absence of degeneracy in one dimension, when a parameter, A, of a potential is varied 
slowly the discrete energy eigenvalue curves, E n { A), cannot cross but they are allowed to come 
quite close and diverge from each other. This phenomena is called avoided crossing of energy levels. 
The parametric evolution of eigenvalues of the generally known one dimensional potentials do not 
display avoided crossings and on the other hand some complicated and analytically unsolvable 
models do exhibit this. Here, we show that this interesting spectral property can be found in 
simple one-dimensional double-wells when width of one of them is varied slowly. 
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In two or more dimensions, evolution of energy eigenvalues as one parameter (A) of the 
system is varied slowly (adiabatically) may show two levels crossing (Fig. 1(a)) or coming 
quite close to each other and then diverging on either side of some special values of A (see Fig. 
1(b)). The former indicates degeneracy (equality of eigenvalues of two linearly independent 
states of one Hamiltonian), the latter is called avoided crossing (AC). In a system, more 
frequent occurrence of AC is debated to be the signature of quantum Chaos [1] and non- 
integrability of the Hamiltonian (non-separability of a Hamiltonian, H(x,y, z..), in various 
dimensions). For instance, the eigenvalues of a particle in a two dimensional rectangular 
box as a function of the ratio of length to breath (L/B) show crossings as its Hamiltonian is 
separable in x and y. But if one small corner of this rectangle is snipped off, the eigenvalues 
show avoided crossings [2] because Hamiltonian is no more separable in x & y or in some 
other special co-ordinates: £,(x,y) & r](x,y). Avoided crossings are known to be responsible 
for surprisingly low energy Landau-Zener [3] transitions. 

Crossings and avoided crossings of eigenvalues can also be understood in terms of time 
independent perturbation theory at the textbook level. When a perturbation is switched 
on, the degeneracy of levels (when A = 0) is lifted and they split into two or more levels (see 
Fig. 1(a)). Similarly, two un-perturbed levels E x and E 2 may split into sub-multiplets. Sub- 
multiplets of different levels may cross (see Fig 1(a)) each other. However, the sub-multiplets 
of one level have the tendency to avoid and diverge from each other as the parameter A is 
varied (see Fig. 1(b)). The Zeeman splitting [3] of different total angular momentum(J) 
states under the magnetic field is one such example. The levels showing avoided crossing 
(Fig. 1(b)) can also called nearly-degenerate(nearly equal) levels and a perturbation widens 
the gap between them. 

On the other hand, the evolution of energy levels of one-dimensional potentials under a 
slow change of the parameter is generally found to be uninteresting wherein all the levels 
increase/decrease monotonically without displaying avoided crossings or crossings (due to 
absence of degeneracy). For instance see Fig. 1(c) for the evolution of first four energy 
eigenvalues of infinitely deep potential of width nX. In this regard, the sub-barrier( below 
the barrier) energy levels of a symmetric double-well [3] potential are interesting which 
display the well known merging of two levels in to one as the height (Vo) or the width ( d ) of 
the in-barrier is increased. For moderate values of Vq and d the sub-barrier levels are closely 
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FIG. 1: Schematic evolution of two eigenvalues depicting (a): crossing of two levels (degeneracy), 
(b): avoided crossing (level repulsion) , (c): evolution of energy levels of one dimensional deep 
infinite well with parameter A: E n = (n/A) 2 . Notice that there is no crossing or avoided crossing 
of two levels. 

lying doublets as in the classical case of Ammonia Molecule [3]. 

Presence of avoided crossings of energy levels in one-dimensional potentials has been ad¬ 
dressed once [5] in extremely complex situations wherein the Schrodinger equation becomes 
Heun equation which can not be solved exactly and analytically. However elegant approxi¬ 
mate treatments have been employed to display avoided crossings of eigenvalues even in one 
dimension. A recent book [6] is a rare avenue where the avoided crossings are well discussed 
and exhibited in piecewise constant one dimensional asymmetric double-well potentials of fi¬ 
nite support in terms of inherent transmission resonances. These potentials of finite support 
represent open systems that can have all: scattering, bound and resonant states. 

Here, we present various analytically solvable asymmetric double well potentials with a 
fixed in-barrier. We show that these closed systems which support only bound states exhibit 
avoided crossing in one dimension when the width of one of the wells is varied slowly. 

In quantum mechanics, we solve the Schrodinger equation: H(x)ip a (x) = E a if; a (x), there 
may or may not be a one to one correspondence between energy eigenvalues E a and ^ a (x). 
For example for free particle case when H = — for one real positive energy E, there 

are two eigenstates e lkx and e ~ lkx , where k = ^ 2 ™' E ■ Similarly in the presence of a potential, 
V(x) for real positive energy continuum there are two states ^(x) and ijj-k(x). 

However in one dimension for bound states of a potential this correspondence is strictly 
one to one as an eigenstate has to satisfy the Dirichlet boundary condition: -0(±oo) = 0. 
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Interestingly, the Schrodinger equation being second order differential equation 

there exist two linearly independent solutions. However, usually only one of them satisfies 
the said boundary condition. Further, it can be proved that even if there are two solutions: 
ip m (x) an d VVi(x) for one fixed discrete energy then they are only trivially different as they 
would be linearly dependent (see Appendix 1). 



FIG. 2: Depiction of various double well potentials, (a) Ml Eq. (2), (b) M2 Eq. (7), (c) M3 Eq. 
(14) and (d) M4 Eq. (19). In the asymmetric harmonic wells the value of the parameter hu on 
the left is more than that of right 

Interestingly, it requires two or more dimensions for degeneracy to occur where energy 
eigenvalues for separable potentials is a function of two positive integers n\ and n 2 . For 
instance in the case of particle in a two dimensional square box of length 7r (2 fi = 1 = h 2 ) } we 
have i.e., E nt n2 = nf+ n|, where E = 260 is four-fold degenerate as 2 2 + 16 2 = 260 = 8 2 + 14 2 
the corresponding eigenstates are if)i(x, y) = A sin 2x sin 16 y and (x, y) = A sin 8x sin 14 y. 
This is called accidental degeneracy. Additional degeneracies arise from the symmetries of 
Hamiltonian. For example, we can have two more eigenstates due to the symmetry of 
Hamiltonian under x -H- y as ij}^{x,y) = A sin 2y sin 16a; and i^ 4 ,(x,y) = A sin 8y sin 14a;. A 
two dimensional Hamiltonian may not be (apparently) separable in x and y then quantum 
numbers ri\ and n 2 loose meaning and eigenvalues are found by semi-classical methods like 
periodic orbit theory [1], 

In more than one dimensions when eigenvalues are studied by varying a parameter 
smoothly and slowly, E num {\) evolve and may become wavy curves rising up and low- 
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ering down having crossings and avoided crossings. The eigenvalues of anisotropic harmonic 
oscillator as a function of anisotropy parameter is a simple example where eigenvalue curves 
are slant straight lines [4], Nilsen’s [4] modified anisotropic Harmonic oscillator model of 
deformed nuclei displays all raising and lowering, crossings and avoided crossings of energy 
levels as a function of the deformation parameter. Eventually, this explains the observed 
nuclear magic numbers. The similar phenomena are also found in quantum dots. 

It is known that the Erst few energy levels of a double well potential are closely lying 
doublets below or around the height of the in-barrier. We bring out a scenario wherein the 
wells of the same depth but different width display avoided crossings when the width of one 
of the wells is varied slowly (adiabatically). Suppose in a symmetric double well potential, 
we have two very close energy eigenvalues, when the width parameter is increased or reduced 
slightly these two levels would diverge from each other to display avoided crossing, what is 
counter intuitive is the occurrence of AC when the asymmetry is large. This is another 
interesting manifestation of Schodinger equation and of quantum mechanics. 

In the following, we End the eigenvalue formulae for four double well potentials depicted 
in Fig. 2. We define u = 2n/Tt 2 and take u = 1 (eFA 02 )' 1 (in Figs. 3,5,6) which corresponds 
to /x ~ 4m e , where m e is the mass of electron in eV, the energy is measured in eV (electron 
volts)and the lengths a, b and c are measured in Angstrom (A 0 ). In the case of Fig. 4, we 
have taken take ji = m e so u = 0.2625 (eFA° 2 ) -1 . 

Ml: Dirac Delta between two rigid walls 
This potential is given as 

{ 0, —a < x < b, 

( 2 ) 

oo x < —a, x > b 

This potential is depicted schematically in Fig. 2(a). When Schrodinger equation has Dirac 
Delta potential at x — 0, the first derivative of the wave function becomes discontinuous at 
x = 0 as [7] 

ip'(x > 0) — ^'(x < 0) = ( 2mV 0 /h 2 )'ip(0 ) (3) 

The appropriate solution of the Schrodinger equation for this double well potential can be 
written as 

{ Asinkix + a), x >—a, 

V ^ (4) 

B sin k(x — b) x < b, 
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satisfying the Dirichlet boundary condition that 'ip(-a) = 0 = 'ip(b ). The continuity and the 
mismatch condition (3) gives 


A sin ka = —B sin kb , 
Bk cos kb — Ak cos ka = (2/i/h 2 )VoAsin ka. 


( 5 ) 


The elimination of A, B gives us the eigenvalue formula for bound states of the potential (2) 
as 

sin/c (a+ 5) = — (2mVo/h 2 ) sin ka sin kb (6) 


When Vo = 0, we get k n (a + b) = rm the energy eigenvalues of infinitely deep well of width 
a + b. When a = b in (6), we get the usual eigenvalue formula for the symmetric case as 
k cot ka = —(mVo/h 2 ) [8]. We take u = l(eVh4° 2 ) -1 , Vo = 10eV, a = 2 A 0 and vary b in 
[0 — 5] A 0 in Eq.(6) to calculate and show the Erst four eigenvalues in Fig. 3(a). Notice ACs 
at b ~ 1, 2,4(A°). The first two ACs are very clear in the enlarged plots 3(b,c). 

In Fig. 4(a), we present the first four eigenvalues E n for Ml for u = 0.2625(eW4° 2 )~ 1 , Vo = 
20eV, a = 5A°. We vary b in [0, 25]A°. Avoided crossings can be seen for E ~ 5,10,15(eF). 
In part (b) we plot E' n = u(a + b) 2 E n to notice the interesting wavy evolution of eigenvalues 
exhibiting ACs clearly. This evolution is similar to the case of Fig. 28.2 and 28.3 of Ref. [6] 
which have been obtained for double wells of finite support and discussed at par with the 
real part of the resonant energies. 

M2: Asymmetric rectangular double well potential 

The simple (symmetric) double well potential is an often-discussed problem in textbooks 
[9,10]. However here we need to discuss the asymmetric double well 


oo, 

x < —a, x > c 


o, 

—a < x < —b , b < x < c, 

( 7 ) 

Vo, 

—b<x<b } 



and obtain the formula for discrete energy eigenvalues. We write the solution of (1) for this 
potential as 


i^(x) = Asink(x + a), x G [—a, —6] 
if}{x) — D sin k(x — c), x E [b, c\. 


( 8 ) 
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For the region —b<x<b we have 


ip{x) = B sinhpx + C cosh pa:, where p = 


\/2m(Vo - E) 


h 2 


We match the solutions and their derivative at x = —b, we get 

A sin k(a — b) = — B sinhpb + C cosh pb 
kA cos k(a — b) = pB cosh pb — pC sinhpb 

Similarly, the matching conditions at x = b give 

B sinh pb + C cosh pb = D sin k(b — c) 
pB cosh pb + pC sinh pb = Dk cos k{b — c) 


(9) 


( 10 ) 


( 11 ) 


We now demand the consistency of the above four Eqs. (10,11) and for their non-trivial 
(non-zero) solutions for A, B,C,D , we get 


— sinhpb cosh pb 0 sin kd\ 

p cosh pb —p sinh pb 0 k cos kd\ 

sinhpb cosh pb sin kd- 2 0 

p cosh pb p sinh pb — k cos kd 2 0 


0 . 


( 12 ) 


Opening the determinant (12) and introducing a — b = d±, c — b = d 2 , d = d\ + d 2 , we get 


y/Eiy o — E) cosh 2 pb sin kd + [E cos kd + V 0 sin kdi sin kd 2 \ sinh 2 bp = 0, (13) 

when E ^ 0, Vo as these cases require separate(special) linear solution ijj(x) = Ax + b [11] 
in the region where E — V (. x ) is zero. 

In Fig. 5(a), we calculate and show the first five eigenvalues of this model (7) for u = 
1 (eW4° 2 ) -1 , Vo = 10 eV,a = 2A°,b = 1A° from Eq. (13). The parameter c is varied 
smoothly in [1, 6]7l 0 . Notice very interesting ACs about the horizontal line E ~ 5.37eE at 
E ~ 2.1, 3.4, 5, 6(eV) which show clearly in the enlarged plot 5(b). 

M3: Dirac Delta in asymmetric Flarmonic well 

This potential can be expressed as 


V(x) = V 0 5(x) 


\pu\x 2 , x < 0, 

|/ic J 2 x 2 ,, x > 0 


(14) 
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(a) 




FIG. 3: (a) Evolution of first four discrete eigenvalues of the model Ml. Vo = 10 eV,a = 2A° 
and b is varied in [0,5]. In (a) for a < 1A° two avoided crossing can be seen clearly. Around the 
symmetric case (a = 2A 0 = b ) there are two avoided crossings. In the enlarged parts (b) and (c) 
the first two ACs appear clearly. Here A = b, \ s = 2A° 

For the harmonic oscillator potential V(x ) = |/xcn 2 x 2 we introduce z = bx, a = s/2jjwjh 
and v = E/(hco) — 1/2 to write (1) as 

' ^ + [u + 1/2 — £ 2 /4]-0(z) = 0, z = bx, a = y/2fiu/h, v = E/(Tuu) — 1/2. (15) 

This equation admits two linearly independent solutions D u (z ), D „(— z) called parabolic 
cylindrical functions. As z oo, D u (z) tends to zero. So in order to satisfy Dirichlet 
boundary condition: -0(±oo) = 0, we choose 

{ AD Ul (—a±x), x < 0, , - 

Zi = aiX,ai = \/2iiUi/h,i = 1,2. (16) 

BD U2 (at 2 x) x > 0 

By matching these solutions and mismatching (3) their derivative at x — 0, we get the 
energy eigenvalue equation as 

+ = (u) 

which simplifies to 

v/ ^ r[ r[-z/ 2 /2] /2] + v/ ^ T[ r[-u!/2} 2] = (18) 

by using the analytic expressions of -0^(0) and D' u ( 0) [12], In Fig. 6(a) we present the 
evolution of first five eigenvalues of the model M3 using (18) for u — 1 (elA4° 2 ) -1 , V 0 = 
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lOeV, Tiuj\ = 2eV by varying hu> 2 in (0, 3]eV. Notice two ACs below Uuj 2 = leV, one above 
E = leV and two ACs around Uu 2 = 2eV (at this value the double well potential (14) is 
symmetric). 

M4: Rectangular barrier in asymmetric harmonic well 

( 


+ a) 2 , 


x < —a, 


V(x) = { 


to, —a < x < a, 

\fi(.A\{x — a) 2 , > a, 

The solution of (1) for (19) can be written as 

' 

FD Ul [ai(—x — a)\, x <—a, 

t/ ’(x ) = A sinh qx + B cosh qx, —a < x < a, 

GD U2 [a 2 (x — a)], x > a, 

Matching the wave function and its derivative at x = —a, we get 

— A sinh qa + B cosh qa = FD Ul (0), 
qA cosh qa — qB sinh qa = —aiFD' (0) 

and at x = a, we get 

A sinh qa + B cosh qa = GD U2 (0), 
qA cosh qa + qB sinh qa = a 2 GD' Uo (0). 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Finding AjB from Eqs.(21,22) equating them and then using the expressions for D v { 0) 
and D' v { 0), we obtain 


r|1 r / [ 2 _ 1/l/ t, )( 21 + ^ r | 1 r f _„ 2 ))( 21 ) 1 ^ h 2 qa = 

\/2/(iVii — E) 


(23) 


e + 2v 4^ rll/2 r[ _^^[^ /2 p /21 ) sinh 2,a, ,= 


r 


Using this Eq. (23), we calculate the evolution of the first five eigenvalues E n ( A) of the 
model M4 for u — 1 (eUA° 2 ) -1 , V 0 = 10eV, fiuj 1 = 2eU and by varying A = Tilu 2 in (0, 3]eV. 
Fig. 6(b) shows straight lines having two ACs below hui 2 = leU, one above E = l(eU) and 
two around the symmetric case hu) 2 = 2(eU). The straight line behaviour of E n ( A) is similar 
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(a) 


(b) 




FIG. 4: First four eigenvalues (E n ( A)) of the model Ml in (a) and effective eigenvalues E' n = 
2+ b) 2 E n /h 2 in (b). We take // = m e ,u = 0.2625 (eVW° 2 ) _1 , Vo = 20eV,a = 5A °, b is varied 
in [0, 25]H°. Here gain A = b, \ s = 5. 



FIG. 5: (a) Evolution of first five discrete eigenvalues of rectangular double well potential (7). We 
take Vo = 10eV, a = 2 A°,b = IA°, and the width of the right well is being varied by changing 
the value of c in [1, 6]. In (a) it appears that the straight line at E ~ 5.37eF is crossing other 
three curves around c ~ 2.1, 3.4, 5,6. But the difference in the types of curves above and below 
this line at negate crossings. In the enlarged figure (b) see the three avoided crossings around 
c = 2.1, 3.4, 5, 6. Here A = c, X s = 2A°. 

to the one obtained in Figs. 1,2 of Ref. [5] where the Schrodinger equation became Heun’s 
second order differential equation. 

Based on our investigations of four double well potentials of various kinds, presented here, 
we conclude that when width (A) of one of the wells is varied slowly the avoided crossings 
of eigenvalues are observed. With the so far understanding of double well potentials, the 
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FIG. 6: First five eigenvalues of the models M3 and M4 when Vo = lOeV in (a) and (b), respectively. 
The parameter twj\ = 2eV for the left well and vary hu 2 in (0,3] for the right well. Notice two 
ACs below huj 2 = leV, one above this and two around Tux 2 = leV. Here A = huj 2 and A s = 2eV 

ACs around the point of symmetry (A = A s ) may be expected. Nevertheless, the ACs much 
below or much above this point are counter intuitive and they invite attention for further 
investigations. 


Appendix 1 


Proposition 1: Let ip m (x), i/j n (x) be two L 2 —integrablc solutions of one dimensional time- 
independent Schrodinger equation satisfying Dirichlet boundary condition: ^(± 00 ) with an 
equal eigenvalue E , then and are linearly dependent. 

Proof: Let the potential V(x) (real or complex) in Schrodinger (2/j = 1 = h 2 ) equation 
gives rise to two solutions ^ m {x) and V ; n(^) with the same energy eigenvalue E, then we 
write 


d 2 ^ m (x) 

dx 2 

d 2 j>n(x 

dx 2 


+ [E - V(x)]ip m (x) = 0, 
- + [E - V{x)]ij> n (x) = 0. 


(A-l) 

(A-2) 


Multiply the first by ijj n (x) and the second by and by subtracting them we get 

d 2 ^ n (x) d 2 i> m {x) df di> n di/j m \ 

^ x) ^~ ~ = C" -*• *T ) = °- (A - 3) 

leading to 


,di> n (x) d^ m (x)\ 

Ym\X) --- i) n {x) --- = C, 


dx 


dx 


(A-4) 
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where C is constant independent of x which can as well be determined at x — ±oo. As the 
eigenstates satisfy ipj(±oo) = 0, we get C — 0. Further we get, implying linear 

independence: (x) = C'Q n (x). 
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